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We determine the detailed thermodynamic behavior of vortices in the 0(2) scalar model in 2D 
and of global monopoles in the 0(3) model in 3D. We construct new numerical techniques, based 
on cluster decomposition algorithms, to analyze the point defect configurations. We find that these 
criteria produce results for the Kosterlitz-Thouless temperature in agreement with a topological 
transition between a polarizable insulator and a conductor, at which free topological charges appear 
in the system. For global monopoles we find no pair unbinding transition. Instead a transition to a 
dense state where pairs are no longer distinguishable occurs at T < Tc, without leading to long range 
disorder. We produce both extensive numerical evidence of this behavior as well as a semi-analytic 
treatment of the partition function for defects. General expectations for N = D > 3 are drawn, 
based on the observed behavior. 
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I. INTRODUCTION 

The many roles played by non-perturbative topologi- 
cal excitations in the dynamics and thermodynamics of 
statistical models and field theories is one of the most 
fascinating and largely unresolved issues of many-body 
systems. 

In the simplest Abelian gauge field theories vortices 
are intimately connected with the existence of a super- 
conducting state in type II materials Q . Excitations car- 
rying topological numbers (instantons, monopoles, vor- 
tices) are also thought to be the best candidates for an 
explanation of confinement in non- Abelian gauge theo- 
ries . In the context of statistical models with global 
symmetries it has long been understood that topological 
excitations can lead to the onset of phase disorder at high 
temperatures Their presence in configurations in 1 
spatial dimensions (ID) in models with local interactions 
down to T = 0+ prohibits, in fact, the establishment of 
long range order at any finite temperature iQ,^. In dis- 
sipative dynamical systems the long range disorder and 
temporal scaling present in the long time limit of phase 
ordering kinetics can also be understood in terms of topo- 
logical excitations 

Several canonical examples illustrate the role of topo- 
logical excitations in bringing about phase transitions 
|^,||J|]. Among them the best known is the Kosterlitz- 
Thouless (KT) transition. At low temperatures the 0(2) 
model in 2D exhibits algebraic long range order, which 
would persist to all temperatures in the absence of topo- 
logical excitations. The advent of disorder at high tem- 
peratures is due to vortex excitations, which appear as 
free charges at the Kosterlitz-Thouless temperature Tkt- 

Recently, due largely to increases in computational 



power and improved methods, many of these systems 
have become available to direct quantitative study. This 
is particularly true of models with global symmetries, 
embodied to a large extent by 0{N) symmetric magnets 
or field theories. 

An interesting question then is: what happens as we 
progressively stray away from cases where topological ex- 
citations are known to dominate critical thermodynamic 
behavior? In this article we investigate this question for 
point defects in 0{N) models in D dimensions, taking as 
a starting point the Kosterlitz-Thouless transition, i.e. 
the case D = N = 2 . The successes of the renormaliza- 
tion group at characterizing critical behavior in D > 2 
suggest that physics in the critical region is dominated 
by perturbative excitations (spin waves). In particular 
as N increases, mean field descriptions become suitable. 
In this limit thermodynamic effects due to topological 
defects are totally unexpected. What then becomes of 
the topological excitations ? Do they disappear from the 
spectrum as likely fluctuations, or do they still occur but 
in a manner that does not lead to long range disorder ? 

The answers to these questions are necessary under- 
pinnings for a general picture of the behavior of topo- 
logical excitations both in equilibrium and as seeds for 
the formation of topological defects upon cooling. The 
current understanding of the formation and evolution of 
topological defects in cosmology and in condensed 

matter requires the presence in the disordered phase of 
fluctuations, which upon cooling can result in long lived 
topological defects. Familiar examples are long vortex 
strings (cosmic strings) or well separated monopole/anti- 
monopole pairs. If these configurations are rare in ther- 
mal equilibrium, above the transition, their abundances 
will be very small and short lived upon cooling. Such 



behavior may have significant phcnomcnological imphca- 
tions and shed new hght on old questions such as the 
monopole problem in cosmology p2| or the planning of 
defect formation experiments in condensed matter sys- 
tems. 

In this paper we study in detail the similarities and 
differences between the statistical behavior of vortices in 
the 0(2) 2D model and of global monopoles in the 'dD 
0(3) model, in thermal equilibrium. For reasons that we 
make clear in Section || this step, between N = D = 2to 
N = D = 3 straddles the boundary between a case where 
topological excitations drive an order-disorder transition 
(the former) and a case where topological excitations may 
be expected to become thermodynamically irrelevant, at 
least for the long distance physics that is characteristic 
of criticality. 

This problem has been investigated in several instances 
in the past, leading to important insights, but a consis- 
tent picture of the thermodynamics of global monopoles 
is yet to emerge. The strongest evidence for an im- 
portant role played by monopoles at criticality in the 
0(3) model in 3D, comes from the study of modified 
partition functions, which include monopole suppression 
terms |13 1^. Lau and Dasgupta |13| showed that the 
introduction of one such term, suppressing all monopole 
fluctuations, results in the disappearance of critical be- 
havior altogether. Later Kamal and Murthy jl4| used 
a different monopole suppression term, which allowed 
monopole anti-monopole pairs with lattice space sepa- 
ration only. They found a new second order transi- 
tion, with exponents different from those of the conven- 
tional 0(3) universality class. Lau and Dasgupta 
also claimed that at the critical temperature of the con- 
ventional 0(3) model the temperature derivative of the 
monopole density dp/dT exhibits a divergence, which 
they argued would signal monopole/anti-monopole pair 
unbinding. This claim was also taken up by Huang, 
Koike and Polonyi who conjectured that the phase 
transition in the 0(3) model would then be driven by 
monopole/anti-monopole separation, in analogy with the 
vortex unbinding that triggers the Kosterlitz-Thouless 
transition in the 0(2) model in 2D. Later the evi- 
dence for a diverging dp/dT disappeared with a high 
precision cluster algorithm study by Holme and Janke 
| p6[ , who showed that dp/dT behaves like the specific 
heat, which does not diverge at Tc- Moreover, Bitar 
and Manousakis searched for unbound monopoles 
by considering phase correlations along closed loops in 
space. They concluded that no such configurations could 
ever be found, implying that the unbinding of monopoles 
plays no role in the critical thermodynamics of the 0(3) 
model in D = 3. 

This body of evidence paints a complex picture of the 
behavior of monopoles at criticality in the 0(3) magnet 
in 31?. It suggests that while monopole degrees of free- 
dom are important in bringing about disorder with in- 



creasing temperature and contribute non-trivially to the 
physics of the iV = 3, I? = 3 universality class, they are 
not essential for the establishment of long range disorder. 
In particular their behavior may not be critical at all at 
Tc. Thus, drawing analogies with the Kosterlitz-Thouless 
transition may give a poor guide to their thermodynam- 
ics. 

The present paper is dedicated to elucidating some of 
these questions, through the detailed comparative study 
of the critical behavior of monopoles and vortices in the 
0(3) Heisenberg magnet and the N = D = 2 XY model, 
respectively. This article is organized as follows. In Sec- 
tion H we review and extend standard free energetic ar- 
guments for point defects of 0{N) scalar models in D 
dimensions. These arguments are both simple and very 
powerful in determining whether topological transitions 
can occur and in elucidating their nature. In Section III 
we characterize the thermodynamic behavior of vortices 
in the Kosterlitz-Thouless transition, by analyzing their 
statistical clustering properties. We find, in agreement 
with the Kosterlitz-Thouless paradigm, that the transi- 
tion proceeds by pair unbinding, which can be observed 
prior to a vortex percolation transition. The latter occurs 
at a slightly higher temperature, approximately where 
the specific heat peaks. Armed with this quantitative 
information of the KT transition and the analytical ar- 
guments of Section ||, we analyze, in Section |^ the sta- 
tistical properties of monopoles in the 0(3) model in 3D. 
We find no unbinding transition as expected on energetic 
grounds alone, but still a percolation transition occurs, 
at a temperature below Tc- We develop methods to de- 
scribe the monopole behavior quantitatively and argue 
that the observed percolation transition can occur with- 
out leading to long range disorder of the order parameter. 
We thus establish the separation between the non-trivial 
monopole thermodynamic behavior and criticality. Fi- 
nally in Section ^ we discuss our results in the larger 
context of 0{N) scalar theories in D dimensions. We ar- 
gue that the criteria of Section II are sufficient to deter- 
mine whether topological defects undergo an unbinding 
transition. 



II. FREE ENERGY CONSIDERATIONS FOR 
0{N) POINT DEFECTS IN D-DIMENSIONS 

We can gain insight into the importance of topologi- 
cal excitations in 0{N) models in D dimensions as vehi- 
cles of disorder by estimating the free energy associated 
with new pair excitations. Later we will specialize to two 
particular cases, those of vortices in 2D and of global 
monopoles in 3D, whose thermodynamics we investigate 
in detail in Section [H . The line of argument used in this 
section follows the original reasoning by Kosterlitz and 
Thouless in motivating the topological transition in 
the 2D XY model, with appropriate generalizations. 



To be definite we consider a general 0(iV)-symmctric 
\(f)^ theory in _D-spatial dimensions. The Haniiltonian is 
written as 



where (f>{x.) is a A^-component real field and \(f>\'^ — 

The 0{N) symmetry of (0) breaks spontaneously at 
low temperatures to 0{N — 1), and the field acquires a 
non-zero expectation value. The degenerate set of min- 
ima lie on a Sn-i sphere. It follows that the homotopy 
group 7rAr_i(S'Ar_i) = Z , the group of integers, which im- 
plies that topological solutions with integer charge exist 
in the spectrum of the theory. In D ^ N dimensions 
these are point defects. The best know cases are the 
kink (or domain wall) in Z3 = = 1, the global vortex 
for I? = iV = 2 and the global monopole (or hedgehog) 
in D = iV = 3. 

These topological defects are classical static solutions, 
i.e. they are (local) energy minima satisfying SH/Scj) — 0. 
In D dimensions point defects arc radially symmetric so- 
lutions. The integer topological charge of these config- 
urations implies a singularity at their origin (r — 0), 
which forces the field amplitude (p{r ^ 0) — > 0. For large 
r it is energetically necessary that the field amplitude 
approaches the minimum of the potential (p(r — > 0) — > 
fo = V- 

These boundary conditions do not guarantee that the 
defect's energy is finite. In fact for _D > 2 the energy 
of topological defects still diverges in the infinite volume 
limit as a consequence of the phase gradient terms, which 
dominate the energy far from the singularity at r = 0. 
These phase gradients lead to an asymptotic form of the 
energy, for large I 



|x|<; 



(2) 



where Uq is the topological charge of the field configura- 
tion. 

The diverging energy of a single defect ioi D > 2 
prohibits it from occurring as a fiuctuation in thermal 
equilibrium in the infinite volume limit. Instead, defects 
can occur in defect/anti-defect multipoles (usually pairs), 
which due to mutual screening can then have finite en- 
ergy, a continuous function of their separation. In the 
case of a pair, the charge separation introduces a natural 
cutoff to Eq. (H) which can be used as an estimate for 
the energy of a pair of size I. This naive expectation may 
be changed for N = D > 3 where the minimal en- 
ergy configuration between two topological defects was 
argued to be one in which the far field is rotated to a 
uniform phase everywhere in space, apart from the re- 
gion between the defects, where energy is concentrated 



and which behaves as a string |20|. Then the interac- 
tion potential between two point topological defects in 
D = N > 3 will be of the same qualitative form as in 
D = N = 3, although the associated string tension will 
differ quantitatively (it is expected to increase with N). 

The simplest interesting example of interacting point 
defects is that of vortices in D = 2. The vortex anti- 
vortex dipole has a field 
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where / is the vector connecting the positive to the neg- 
ative charge in the pair. As a consequence of (^) a point 
vortex far away from the dipole interacts with it via a 
potential inversely proportional to distance. Two well 
separated pairs then interact with a potential 



(4) 



where li and I2 are the separation vectors within each of 
the pairs, r is the vector connecting the center of the two 
pairs and f — r/r. Thus a dilute gas of weakly interacting 
vortex pairs can exist at low temperatures. 

Global monopoles have stronger linearly confining po- 
tentials. Their dipoles therefore behave as 
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Although finite the potential of this dipole is still sub- 
stantial. Two well separated pairs of monopoles then 
interact via 



h.h-{f.h){f.h) ^^^P^ 



(6) 



Since the interacting potential decreases inversely with 
separation we see that monopole pairs can similarly ex- 
ist in a dilute, weakly interacting state, but also that 
their mutual interactions are stronger than between vor- 
tex pairs. As we will see below this characteristic af- 
fects the thermodynamics of monopoles relative to vor- 
tices considerably. We remark in passing that the leading 
interaction between monopole pairs, apart from polar- 
ization inner products, behaves in 3-D as the Coulomb 
potential. The Coulomb gas in 3-D has a transition, as- 
sociated with the familiar process of ionization, that is a 
smooth analytical crossover. In terms of monopole pairs 
this transition would occur between a phase of free pairs 
gas and another where monopole pairs become bound 



to form clusters. If it occurs, we may then expect that 
this monopole pair transition wiU not lead to critical (i.e. 
non-analytic) behavior. 

To estimate the free energy of a pair of defects we must 
finally estimate its entropy S = log{n), where fl is the 
number of states of the pair, is proportional to the 
surface of the D — 1 sphere of radius I i.e. the number 
of configurations a pair can take when rotated around its 
center of mass. 

Then the single pair free energy in arbitrary dimension 
D is 

Fz3=2(0- Ec + ao log(0 -Tho log(0, 

Fz5>2(0- + an /^"^ -T{D-l)bD log(/), (7) 

where Ec accounts for the total core energy of the two 
defects in the pair and ajj and b d geometrical dimension- 
less constants dependent on space dimension. 

We emphasize that in these considerations we ne- 
glected the eff'ects of other defect pairs. Qualitatively 
these will reduce the free energy of the new pair rela- 
tive to the above estimates, as they will tend to orient 
themselves in order to (partially) screen the new charges. 
Thus (0) should be thought of as an overestimate. 

The free energy of the pair gives us a qualitative mea- 
sure of its probability in equilibrium P{1) cx e^^^'-'/"^. 
We will explore this relationship further in Sec. III. For 
now we note that for D = 2 both the energy and entropy 
terms behave logarithmically with I and the overall sign 
of the free energy for large pairs depends on the tem- 
perature, as noticed by Kosteiiitz and Thouless ||l^. In 
the low-T regime F grows with charge separation leading 
to suppression of large pairs. For high temperatures the 
entropy term is dominant and large pairs have negative 
free energy. This suggests the existence of a high-T phase 
characterized by the unbinding of defect/anti-defect pairs 
and the production of free charges. This is the essential 
idea behind the Kosterlitz-Thouless mechanism for the 
0(2) 2D transition. 

For higher dimensions D > 2 the energy term domi- 
nates at large I for all temperatures and large pairs always 
remain exponentially suppressed. This simple argument 
suggests that there is no unbinding topological transi- 
tion for I? > 2 in global models and that defects remain 
tightly bound for all T. In order to destroy this picture it 
is necessary that the behavior of the bare energy and/or 
number of configurations with I would change due to in- 
teractions with other defect multipoles. Such behavior is 
not seen in the Kosterlitz-Thouless transition where the 
effect of other pairs softens the field modulus ipo (the su- 
perfluid density or spin wave rigidity) , leading to a renor- 
malized value of the transition temperature but not to a 
different kind of transition. In Eq. (|^) we have taken the 
pair energy to be determined by a simple cutoff of the 
single charge total energy (^. Our conclusions would 
remain unchanged if as suggested in [n9| for D > 2 the 



interaction becomes linear as the field's interacting core 
collapses into a string connecting the two charges in the 
pair. The entropy term would still be dominated by the 
interaction energy in the same way as for D = 3. 

Thus we have reached the expectation that no unbind- 
ing topological transition should ever occur for N = D > 
2. We put this expectation to the test in Sec. Ill, where 
we study comparatively the thermodynamic behavior of 
both vortices in 2D and global monopoles in 3D. 

We devote the remainder of this section to a few addi- 
tional remarks about the applicability of the free energy 
considerations of defect pairs to more general circum- 
stances. An interesting limit is that of systems that re- 
main disordered due to topological configurations down 
to T = 0. In the class of models of (|l|) only the case of 
the = 1 A^"* theory (or the Ising model) in ID has 
this property, due to the presence of kinks (or domain 
walls) . In gauge-t-Higgs field theories the physical proper- 
ties of topological solutions change radically because the 
phase gradients, which dominate the energetics of global 
defects, become pure gauge transformations and carry 
therefore no energy. This property is a direct result of 
the Higgs mechanism. The phase gradients correspond to 
Nambu-Goldstone modes, each associated with a genera- 
tor for the remaining unbroken symmetries. In the Higgs 
mechanism these massless modes are 'eaten' by the gauge 
field, which in turn acquires a mass for its longitudinal 
projection. Thus the total energy of gauge defects is con- 
centrated in their cores, and falls off exponentially with 
distance. Then gauge topological charges interact via a 
short range potential, in contrast to global defects. This 
interaction energy can typically be neglected or treated 
effectively, as a change of the core energy, in our thermo- 
dynamic estimates. 

Either way by repeating the free energy argument with 
core terms only we see that the entropy contribution is 
dominant for all temperatures, for large separations in 
D > 2. This implies that for large enough separation 
defect/anti-defect pairs are always likely fluctuations, 
and suffer no Boltzmann suppression down to T = 0^. 
This 'condensation' of free topological excitations can 
explain striking properties of non-Abelian gauge theo- 
ries pl| . The thermodynamic spectrum of these mod- 
els should then be characterized by the existence of a 
dilute gas of free defects at low temperatures. The nu- 
merical verification of this expectation is presently the 
subject of intense research. Hints of this behavior have 
recently been found numerically, see e.g. for the case 
of a SU{2) lattice-gauge theory in 31?, which possesses 
't-Hooft-Polyakov magnetic monopoles - the direct ana- 
logues of the global monopoles in (|^) for A^ = Z? = 3. 



III. CHARGE CLUSTERING AND THE 
KOSTERLITZ-THOULESS TRANSITION 

As outlined in Section |l| the behavior of a system 
of charges with logarithmic interactions in 2D (the 2D 
Coulomb gas) is well known, underpinning the topolog- 
ical transition in the 0(2) 2D model. In this section 
we develop diagnostics that allow us to measure, in nu- 
merical studies of the equilibrium partition function, the 
critical properties of topological point charges. Later we 
will use the Kosterlitz-Thouless behavior as the bench- 
mark for a charge unbinding transition of iV = _D = 3 
monopoles. 

Rather than using a discretized two dimensional ver- 
sion of the Hamiltonian in Eq.(|]) we chose to study the 
2D XY model which belongs to the same universality 
class. This choice has the advantage that the XY model 
thermodynamics has been extensively studied, both an- 
alytically and via large scale numerical simulation. Con- 
sequently its critical properties are well known quanti- 
tatively, including dimensionful quantities such as the 
Kosterlitz-Thouless transition temperature Tkt and the 
temperature at which the specific heat peaks Tcv- 
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FIG. 1. The plaquette density of vortex pairs for the 2D 
XY model. Error bars denote standard deviation over 4000 
independent field realizations. Both the Kosterlitz-Thouless 
temperature Tkt and Tcv, at which the specific heat peaks, 
are shown. No signs of critical behavior are apparent in the 
total vortex density. 

The XY model consists of a set of two-dimensional 
unit-length spins with nearest neighbor interaction. The 
Hamiltonian is given by: 



H[{s}] 



(ij) 



(8) 



where the sum is over all pairs of nearest neighbor sites 
and we take J — I. 

All quantities below were obtained via standard 
Monte-Carlo generation of large ensembles on a lattice 
of size 128^. For each temperature we obtained a set 
of 2000 independent configurations, from which we mea- 
sured global properties of the vortex population. Lo- 
cal quantities not involving use of time-consuming clus- 
ter algorithms were averaged over larger ensembles. The 
vortex content of each field realization is determined by 
identifying integer spin windings around the lattice pla- 
quettes. The values for the two characteristic tempera- 
tures, Tkt = 0.89 and Tcv = 1-03 were obtained from 
the literature (see e.g. |2^j2^) and are confirmed below. 
In particular we checked that the specific heat peaks at 
Tcv within statistical error. 

Fig. 1^ shows the temperature dependent density of vor- 
tex pairs pvv (T) defined as the fraction of lattice plaque- 
ttes occupied by a positive charge. Note that although 
the total vortex density pvv(T) increases with tempera- 
ture it does not show a clear change of behavior at either 
Tkt or Tcv- This is not surprising since the system does 
not undergo a 2"'^-order phase transition, and the criti- 
cal singularities are much weaker in nature. In particular 
the properties of a few large vortex pairs, crucial for the 
onset of phase disorder, are masked in pvv(T) by the ex- 
istence of many more small pairs. In order to see signs of 
the unbinding we must study the properties of the vortex 
population in greater detail. 

To achieve this we must deal with the ambiguity in- 
volved in grouping vortices in pairs. To overcome this 
problem in the most general way possible we choose to 
group the vortices in each field realization into clusters, 
defined in terms of an adjustable length parameter Zci- 
Vortices or anti-vortices - we do not distinguish between 
the two - separated by less than Zci are collected in the 
same cluster. Thus, each cluster consists of the set of 
vortices and anti- vortices that lie within a distance lc\ of 
at least another element of the cluster. 

The cluster decomposition is achieved efficiently by 
applying a generalization of the Hoshen-Kopelman algo- 
rithm [ p5| , developed originally for studies of percolation. 
Id is successively increased, starting from the lattice spac- 
ing, the smallest length scale in the problem. For each 
Ici we measure a set of cluster properties. In particular 
the topological charge properties of clusters are ideal di- 
agnostics in the search for signs of a charge unbinding 
phase transition. 

As a consequence of our choice of periodic boundary 
conditions the sum over the charge of all clusters in the 
volume is always zero. To quantify the typical charge of 
a cluster, we define, for each choice of Id, a mean cluster 
charge Qd- Qd results from adding up the absolute value 
of the charge Q of each cluster in a given realization and 
dividing the result by the total number of clusters, i.e. 
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(9) 
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Fig. H shows Qci(^ci) for two different values of the tem- 
perature around Tkt- 



0.06 

0.04 

0.02 

J 
0.06 

0.04 

o.os 



T = 0.87 



_i — ° ° I ? " T "? " T " f" •y - 1 - T- -T - 



T = 0.95 







10 
1„, 



15 



20 



FIG. 2. Cluster charge Qci vs. clustering length l^i for a 
temperature slightly below (top plot) and above (bottom plot) 
the Kosterlitz-Thouless temperature, in the 2D XY model. 
Below the transition the mean charge decreases monotonically 
with Ici- The clustering length Id for which Qd ~ increases 
as T — > Tj^rp. For T > Tkt, Qci peaks at an intermediate 
value of lc\ before decaying to zero. Standard deviation error 
bars of order of the results not shown for clarity. 

The behavior of Qci at low temperatures can be under- 
stood in terms of the properties of a dilute distribution 
of vortex/anti-vortex pairs. A vortex and an anti-vortex 
will be in the same cluster if their separation is smaller 
than the chosen ^ci- If the distance between different pairs 
is large, all pairs with separation smaller than Id will be 
in neutral clusters. In this case the charged clusters will 
consist only of single charges from pairs with separation 
larger than ^ci- With increasing clustering length, Qd 
should then decrease, reaching zero when ^ci becomes of 
the order of the largest pair in the sample. At this length 
all clusters become neutral. We then define Lnoutrai by 
Qci (-^neutral) = 0, a mcasure of the size of the largest 
pairs at a given temperature. 

The top plot in Fig. || illustrates the typical behav- 
ior of Qci just below Tkt- Qci(^ci) has a long tail for 
large values of the clustering length (up to Id — 20), 
signalling the presence of large pairs. For lower temper- 
atures, Qci decays faster. In Fig.^ we plot the variation 
of ijicutrai with T. Up to the Kosterlitz-Thouless tem- 
perature incutrai iucreases as pairs with larger and larger 
separations are produced. 




FIG. 3. The smallest clustering length Lncutrai, at which 
all vortex clusters are neutral, vs. T, in the 2D XY model. 
For low densities Lncutrai corresponds to the separation of the 
largest pair in the (finite) ensemble. 

According to the standard picture of the KT transition, 
above Tkt free charges appear in the system. Their pres- 
ence affects Qd{lc\) because a population of free vortices 
changes the charge of otherwise neutral clusters. Thus 
we can no longer expect Qci('ci) to decay monotonically. 
In fact we observe that above the transition Qd displays 
a peak at a finite value of the clustering length, which we 
define as Tpoak- With increasing temperature the value 
of Tpcak decreases and the height of the peak increases. 
This is the result of a higher vortex density and, among 
them, more free charges. The increase in the density 
of free vortices also reduces the mean distance between 
them, and moves the peak to lower Id- This behavior has 
the important characteristic that Tpoak diverges as Tkt 
is approached from above. Higher free charge densities 
also lead to a decrease in Tncutiai above Tkt- 

The behavior of both Tnoutrai and Tpcak is shown in 
Fig. H and Fig. ^ respectively. Both quantities show clear 
signs of critical behavior at Tkt 7 although the behavior 
of Lncutrai IS plagued by higher statistical uncertainty. 

Note that pair unbinding is not the only way a set 
of point charges may display a change of properties. In 
general as temperature is increased two concurrent effects 
take place. The first is that pairs with larger separation 
and higher interaction energy are nucleated. The sec- 
ond is the production of more pairs at small separations. 
Depending on the interplay between these two trends a 
situation may be reached when the distance between dif- 
ferent pairs is of the same order as the separation within 
each pair. In this case the system becomes dense (it per- 
colates) and pairs become indistinguishable. 




FIG. 4. Length for cluster charge peak in the 2D XY 
modeL For T > Tkt, the length at which the cluster charge 
function peaks decreases with T due to production of higher 
densities of bound and unbound pairs. Below Tkt the mean 
charge decreases with l^i; we define I/poak = oo in this case. 

In order to determine the temperature at which pair 
percolation occurs in the 2D XY model we measured, 
for each configuration, the value of the clustering length 
at which all clusters become neutral^ (Imax) as well as 
the minimum /ci, for which all vortices in the sample fall 
within the same cluster, -Lpcrc- 

For dense vortex systems, where there is no positive 
correlation between vortices and anti- vortices, (/max) — 
Lpcrc- This limit must be reached at high temperature 
as is indeed shown in Fig. ^. Moreover we see that for 
T < Tcv the length (/max) < ^porc- This includes the 
vicinity of Tkt, where pairs remain dilute enough that 
they can be identified. The density threshold where the 
vortex system becomes dense is T ~ Tcv, but a precise 
identification would demand careful finite volume scaling. 
In any case we also see that the approach to a dense state 
occurs seemingly continuously, without any clear signs of 
critical behavior. In this sense it may not be possible to 
associate it with a particular value of T. 



* (/max) differs from Lneutrai in the sense that (/max) is a 
thermal average of the size of the largest pair in each sample, 
whereas Tneutrai is the size of the largest pair in all samples 
in our ensemble. In this sense the peak in (/max) reflects a 
maximal production of large pairs where the peak in Lncutrai 
corresponds to the production a single very large (presumably 
infinite in the infinite volume limit) pair. 




FIG. 5. Mean clustering length for neutral clusters (/max) 
in the 2D XY model, error bars corresponding to standard 
deviation over 2000 field samples. The dashed curve shows 
the clustering length that determines a single cluster - stan- 
dard deviation of the same order as for (/max), not shown 
for clarity. When the two curves meet near Tcv the system 
becomes dense. 

We have now used cluster decomposition methods ap- 
plied to the vortex population to characterize its critical 
properties. As the temperature is increased we see that 
free vortices first appear at Tkt and are maximally pro- 
duced approximately at Tcv, where the vortex system 
becomes dense and the concept of a vortex pair ceases to 
be meaningful. 



IV. MONOPOLES IN THE 3D 0(3) MODEL 
A. Thermodynamics of the model 

In this section we apply the tools developed in the con- 
text of the Kosterlitz-Thouless transition to a 3D scalar 
field theory with 0(3) symmetry. Our analysis will be 
based on a discretized version of a theory, Eq.(0). 
We start by establishing its thermodynamic properties. 
In particular we will be interested in determining the 
value of the critical temperature Tc at which the model 
displays a second order phase transition. 

In order to generate a Boltzmann distributed ensem- 
ble of field configurations we have evolved a second-order 
in time Langevin field equation (see ||2^ for more de- 
tails). For our purposes this is equivalent to using a 
Monte-Carlo or cluster algorithm. The advantage of the 
Langevin equation is that it can be easily generalized 
for time-dependent systems. In a future publication the 
equilibrium states generated in this way will be taken as 



initial conditions for real time out-of-equilibrium studies 
&■ 

The procedure is as follows. We evolve the three com- 
ponent real scalar field in time with the equation of mo- 
tion 



id'. 



3 = 1:3 



+ +?74=r„ (10) 



where i S {1,2,3}. We discretize this scheme using 
a staggered-leapfrog method with time-step 5t — 0.04. 
The random force Ti{x, t) is a Gaussian distributed field 
with temperature T as determined by the fluctuation- 
dissipation theorem. Ti is characterized by 

{T,{x)) = 0, {T,{x)T,{x')) - |?<5,,(5(x - x'). (11) 

The value of the dissipation coefficient r] does not influ- 
ence the equilibrium results and in our simulations we 
chose ?7 = 1.0 to ensure rapid convergence. The lattice 
spacing was set to 5x — 0.5 and the model parameters 
chosen to be w? = 1.0, A — 1.0. We used a computa- 
tional domain with L = 100 points per linear dimension. 
Both local and global observables were measured over at 
least 200 independent field realizations. 

As an order parameter we have used the norm of the 
spatially averaged field (|'/>y|), defined as: 



(l<^y|> = 



which is analogous to the magnetization in spin models. 




(T-Te)/T, 

FIG. 6. Order parameter for the 3-D 0(3) model and corre- 
sponding power-law fit in the critical region. The temperature 
has been rescaled to (T — Tc)/Tc setting Tc — 0.41 as deter- 
mined from the fit. Error bars denote standard deviation over 
an ensemble of 200 independent field realizations. 



Fig.^ shows the temperature dependence of (|(/>y |) . For 
T > Tc, vanishes. Near but below Tc, the order 

parameter displays universal critical power law behavior 



{\^v{T)\)=B{ 



Tc-T. 



Tc 



r, /3>o. 



(13) 



which is the analog of the magnetization density in spin 
models. Here (3 is the universal critical exponent as- 
sociated with the behavior of the magnetization below 
Tc and is not to be confused to the inverse temperature 
elsewhere. By fitting the numerical values for to 
Eq. ( |l3| ) we are able to measure the critical temperature 
obtaining Tc — 0.41. This sets a reference point, the most 
important scale in the system. We also compute the crit- 
ical exponent f3 = 0.36. This is in good agreement with 
both recent theoretical and large scale Monte-Carlo esti- 
mates for the critical exponent which give f3 = 0.366(2) 
and P = 0.3685(11) respectively (see e.g. and refer- 
ences therein) and provides a check on the accuracy of 
our numerical setup. 



B. Monopole Statistics 

We are now ready to analyze the properties of the equi- 
librium global monopole population. For each field real- 
ization monopoles and anti-monopoles are identified by 
measuring the three dimensional winding number of the 
field around each cubic lattice cell. The algorithm used 
is based on a higher dimensional generalization of the 
geodesic rule traditionally used for identifying strings in 
0(2) theories Details of this procedure are given in 
Appendix H. 
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FIG. 7. The mean plaquette density of monopole pairs pmm 
vs. reduced temperature for the 3D 0(3) model - error bars 
as in Fig. n. 




Fig. 1^ shows the density of monopole/anti-monopole 
pairs (defined as the total positive charge in the compu- 
tational domain divided by the number of lattice sites) 
versus temperature. The total monopole density Pmni 
increases smoothly with T and its derivative peaks at 
the critical point signalling the second-order phase tran- 
sition. Above Tc the rate of increase diminishes and 
the total pair density converges slowly to approximately 
Pmm — > 0.17 as r — !■ cx) (not shown). 

Fig. H shows a log-linear version of Fig. illustrating 
how over nearly the whole temperature range below Tc, 
Pmm is well fit by an exponential Pmm — Ae~^°^'^. Only 
very near the critical point does the fit fail to follow the 
density curve accurately. This behavior suggests that the 
increase in the total monopole density is dominated by 
the creation of large numbers of minimum size pairs, each 
with typical energy Eq ~ 2.0. 
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FIG. 8. Fit for the mean density of monopoles to 
^exp(— i?o/r), error-bars as above. In the low temperature 
regime we obtain A = 7.0± 0.6 and Eq = 1.98 ± 0.02, fitting 
the first 23 data points. Below Tc the density displays near 
exponential behavior down to T = 0. 

An understanding of this behavior can be obtained by 
evaluating the partition function for monopole pairs un- 
der certain simplifications. If we assume pairs are inde- 
pendent i.e. we neglect pair/pair interactions and volume 
exclusion effects, the partition function for a pair is 



-Ep/T 



(14) 



where the sum is taken over all single pair internal con- 
figurations, i.e. it excludes translational modes. We take 
the pair energy to be of the form 



where I is the monopole anti-monopole separation in 
units of lattice spacing. In the continuum the simplest 
way to compute Z{T) would be to use the approximate 
expression for the free energy Eq.(|^). While this should 
be valid for large values of the pair size, such evalua- 
tion of the number of states breaks down relative to that 
on the lattice, especially when I becomes of the order of 
the lattice spacing. Since we expect the monopole pop- 
ulation to be dominated by small pairs, the continuum 
approximation would be a significant source of error. To 
circumvent this problem we calculate the partition func- 
tion by evaluating numerically the sum in ( |l^ over all 
possible lattice configurations of a pair with fixed center, 
for a given choice of Ef. and a. In this way the number 
of pair configurations on the lattice, and therefore the 
entropy, are calculated exactly. 
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FIG. 9. Numerical data for the total pair density and cor- 
responding fits near Tc, error-bars as in Fig. ^. The en- 
ergy varying theoretical curve (solid line) was fitted to the 
18 lowest temperature points (not shown in the plot) up to 
(T - Tc)lTc = -0.25. The simple exponential fit (dashed 
line) fails to follow the numerical curve up to Tc indicating 
relevant production of higher- energy, larger separation pairs 
in the critical region. The varying energy fit on the contrary 
matches the data well up to the critical point. 

Monopole thermodynamic averages can be easily eval- 
uated from Z{T). The total pair density per site is given 
by: 



Pmni(r) 



Z{T) ~ 1 



(16) 



E.„ ^ Ec + (jI, 



(15) 



By allowing Ec and a to take arbitrary values, we fit 
the predicted pair density to the numerical data. A x 
squared fit leads to estimated values of the interaction 
parameters of Ec — 0.26 and a ~ 1.71. 



In Fig. 1^ we compare the monopole density Pmm ob- 
tained in this way to the numerical data and to the sim- 
ple exponential fit discussed before. Clearly taking into 
account pairs with variable separation improves the esti- 
mate leading to precise results up to Tc- In spite of this 
improvement it remains true that the monopole thermo- 
dynamics is always dominated by very small pairs. Us- 
ing the fit results for Ec and a we obtain the follow- 
ing energies for the smallest pairs allowed on the lattice: 
E{1) = 1.97, E{^/2) = 2.68 and E{^/i) = 3.22. Con- 
sidering the contribution of the first three terms in the 
partition function already leads to a very reasonable ap- 
proximation to Z(T): 

Z{T) = 1 + 6 + 12 + 8 e-^-^^'^ . (17) 

This gives /Omffl(Tc) = 0.07, in good agreement with the 
measured value (the integer pre-factors are the number 
of different lattice configurations for a pair at these sep- 
arations). 

The value measured for a is considerably lower than 
the one obtained from the single monopole classical esti- 
mate in Section ^ Evaluating Eq. (^ exactly we have 
for the single pair energy Ep = Ec + Atttit? / X x Z, which 
leads to (7 = 27r in lattice spacing units {Sx = 0.5). The 
difference between this classical value of cr and its value 
inferred from fitting the thermodynamic monopole den- 
sity is a consequence of strong medium dressing effects, 
resulting both from the influence of other monopole pairs 
and from interactions with the spin wave degrees of free- 
dom. 

In any case the unequivocal exponential behavior of the 
total density below the critical point implies that cr < 2. 
The quality of the fit using cr = 1.7 and its success in 
predicting other features (see below) of the data suggests 
that one should not place exaggerated confidence in the 
classical single monopole result. 

In a previous publication ]3C| ] two of us predicted 
the value of defect density at criticality for 0{N) the- 
ories. This calculation assumes that fluctuations at Tc 
are Gaussian, with their scale invariant connected 2- 
point function characterized by the universal critical ex- 
ponent 77, the anomalous dimension. Using the value of 
Pmm — 0-17 at infinite temperature as a normalization 
(see for details) leads to a predicted /Oinm(7c) — 0.07 
for 0(3) in 3-D, in good agreement with the present nu- 
merical measurements. 

A similar calculation of the temperature dependence of 
the vortex pair density can be done for the 2D XY model. 
As in the monopole case, the low temperature Monte 
Carlo data is reasonably well fitted by an exponential. 
Assuming a pair energy of the form Ep = Ec + a log(/) 
and calculating the partition function as before, the pre- 
diction can be improved leading to good results up to 
Tcv (see Fig. ^0|). The single exponential fit for low-T 
gives Eq ~ 7.2 which compares well with a previous result 



of 7.5, measured by Gupta and Baillie The differ- 

ence is probably due to our fit being based on low-T data 
points only. In the same article Gupta and Baillie 
also obtained a different exponential fit in the tempera- 
ture region between Tkt and Tcv with a higher value for 
Eq. Using the form (|l5| ) we are able to fit both temper- 
ature regimes obtaining the x squared results Eq ~ 6.7 
and cr ~ 2.9. 
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FIG. 10. Same as Fig.^for the 2D model. The theoret- 
ical curves were obtained by fitting the low-temperature data 
up to T = 0.75. Whereas the simple exponential (dashed line) 
fails to follow the observed density above Tkt- The curve ob- 
tained from the partition function, which includes pairs with 
all separations, matches the data well up to Tcv- 

Following in the footsteps of the 2T) charge cluster 
analysis we now turn to the properties of monopoles in 
the 0(3) ZD theory. Fig. |l| shows (Zmax) and Tperc in 
terms of the reduced temperature in the critical region. 
The monopole ensemble becomes dense when the two 
length scales are comparable, (Imax) — iperc and it is no 
longer possible to identify isolated pairs. This happens 
at [T — Tc)/Tc ~ —0.25, the temperature at which (Zmax) 
peaks, well below Tc- This behavior stands in striking 
contrast to that of vortices in the 2D case (see Fig]|), 
where the vortex gas percolated only in the exponentially 
disordered phase at T = Tcv > Tkt- 

It is important to realize that this behavior of 
monopoles is not in contradiction with maintaining long 
range order up to Tc. The system of monopole anti- 
monopole pairs can become dense without disordering 
the field over large distances. This can be understood by 
considering a domain with radius much larger than the 
maximal pair size and is essentially the earlier result of 
Bitar and Manousakis |lj] . Since the total field winding 
in the domain's surface is given by the total charge in its 



interior, its value will be zero. That is, finite pairs will 
not affect the long range behavior of the field on scales 
larger than their size, only unbound charges can lead to 
the break down of long range order. 
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FIG. 11. Mean clustering length for neutral clusters (er- 
ror-bars correspond to standard deviation over 200 field sam- 
ples) and percolation length (dashed curve - standard devia- 
tion error bars of the same order as for (/max), not shown for 
clarity) for the 3-D 0(3) model. The system becomes dense 
well below the critical temperature when the two quantities 
become comparable. 

The proliferation of small pairs in the three dimen- 
sional case is made possible by the fact that the core 
energy of the monopole is small when compared with the 
interaction potential. This leads to production of large 
densities of small pairs, while large pairs remain strongly 
suppressed by the fast growing linear interaction term. 
As a consequence the system percolates mostly due to 
a high density population of small pairs at a tempera- 
ture where large pair configurations have a exponentially 
negligible contribution to the thermodynamics. 

To investigate this behavior we can use our approxi- 
mate partition function to calculate N{1), the density of 
pairs of size I as 



N{1) = Z-i(T)n(Oe-^''(''/^, 



(18) 



where n{l) is the total number of configurations for a 
pair of size I. This expression can be readily evaluated 
numerically. We then use N{1) to estimate the size of 
the largest pair in a computational domain by finding 
the value for I such that N{1) x ~ 1., where L is 
the linear size of the computational volume. That is, we 
demand that in each computational volume there should 
be on average one pair of maximal size. This length scale 



corresponds to (/max)- The value of (/max) estimated in 
this way is plotted against temperature for both models 
in Fig. |l|. 

In order to calculate the percolation temperature in 
this approximation we must estimate L 



pcrc ■ 



This can 

be done by assuming that the typical distance between 
pairs is of order of 1/ p*-^/^' , where p denotes either pvv or 
Pmin depending on the dimension. This implies that we 
will have percolation when ~ l/p(^/-^\ Fig. 1^ shows the 
temperature dependence of ipcrc, estimated as l/p^'^/^\ 
using the parameters from the previous density fits. The 
point where the two curves meet defines the percolation 
temperature where the change ensemble becomes dense. 
The values obtained agree reasonably well with the data. 
For the 3D 0(3) case we find the percolation temperature 
to be (T — Tc)/Tc ~ —0.2 with pairs of maximal size 4, 
compared to (T — Tc)/Tc ~ —0.25 with largest length of 
around 5.5 from the numerical results. For the XY model 
we obtain T ~ 0.96 slightly below Tcv but still clearly 
in the exponentially disordered phase and / ~ 10, which 
coincides with the numerical result for the mean size of 
the largest pair per box at the percolation temperature. 
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FIG. 12. Theoretical prediction for the temperature depen- 
dent mean maximal pair size (/max) and pair separation Lperc, 
for the 2D XY model (left) and the 3D 0(3) model (right). 
At the temperature where the two curves cross, the system 
becomes dense and it is no longer possible to distinguish in- 
dividual pairs. The step-like appearance of (/max) in the right 
plot is due to lattice discretization effects on /, which are ap- 
parent due to the small size of the monopole/anti-monopole 
pairs. 



V. CONCLUSIONS 

We started this paper by using single defect pair free 
energetic arguments to predict that no pair unbinding 
transition should occur for point defects in 0{N) models 
with global symmetries, for N = D > 2. We then veri- 
fied this prediction by comparatively studying the defect 
thermodynamics of vortices and global monopoles oc- 
curring in typical configurations drawn numerically from 
canonical ensembles of the 0(2) model in 2D and an 0(3) 
field theory in 3D. 

We measured in great detail the behavior of vortices 
across the Kosterlitz-Thouless transition and found it in 
agreement with the expectations of the theory, which 
predicts the appearance of free vortices at Tkt- In the 
0(3) model in 3D isolated single monopoles never oc- 
cur. Instead the monopole ensemble transits from a di- 
lute pair phase at low temperatures, to a dense monopole 
gas, through the nucleation of large numbers of tightly 
bound pairs below Tc. We have seen that this behav- 
ior is consistent with qualitative expectations based on 
the single monopole pair energetics, which is character- 
ized by a light core and a linearly confining interaction 
potential. Our quantitative treatment shows, however, 
that the bare parameters in the monopole potential, the 
monopole core energy Ec and string tension tr, are sig- 
nificantly renormalized by thermal medium effects. 

The thermodynamic behavior of global monopoles is 
consistent with both a (non-critical) topological charge 
conductor to insulator transition and the absence of long 
range phase disorder at temperatures below Tc- Global 
monopoles in 31? do not behave like Coulomb charges, 
but rather more like static quarks, since they interact 
via a linearly confining potential. The conductor phase 
in this case cannot be reached by the nucleation of truly 
free isolated charges, as in the Kosterlitz-Thouless transi- 
tion. It is instead the result from the fact that in a dense 
monopole gas charges can move freely, by hopping from a 
nearby anti-charge to another. Because local charge neu- 
trality persists to all finite temperatures the dense 
monopole gas does not lead to long range disorder of 
the field phases (spins). Long range phase fluctuations 
must therefore arise from different degrees of freedom, 
the spin waves. This is indeed the renornialization group 
picture of the transition, which correctly predicts all uni- 
versal critical behavior. From this perspective we con- 
clude that monopole excitations, although contributing 
to short range disorder in the 0(3) model in 3-D, are 
incapable by themselves, because of their enslaving en- 
ergetics, of producing long range disorder at any finite 
temperature. This is not to say that their disordering in- 
fluence can be completely neglected. Their effects must 
generate a short distance nontrivial quantitative contri- 
bution to critical exponents, thus making the universality 
class in their absence different . 



In this way we must conclude that there is a non-trivial 
interplay between monopoles and spin waves in the 0(3) 
model and that the detailed nature of critical behavior is 
changed if either are suppressed. 

The generalization of these results to arbitrary N = 
D > 3 is immediate. The potential between defects be- 
comes steeper and steeper as a function of I, and topo- 
logical excitations, just like monopoles in 3-D, will never 
unbind. In this way topological excitations do indeed be- 
come more and more irrelevant, as N increases, for the 
physics of large length scales that characterizes thermo- 
dynamics in the critical domain of 0{N) models. 

We also see from this perspective that, upon cooling 
the system, most topological fluctuations will annihilate 
with a nearby anti-defect, leading to small and quickly 
disappearing populations of topological defects. From 
this perspective global topological monopoles formed at 
a cosmological phase transition present no real danger of 
creating a 'monopole problem'. 

We conclude by invoking a complementary view of the 
phase transition in 0{N) models. There is clear ev- 
idence (and mathematical proofs in certain particular 
cases) that criticality in 0{N) scalar models is equiva- 
lent to the percolation of so-called Wolff spin clusters . 
Wolff clusters are built by forming bonds among adjacent 
spins according to a temperature dependent probability. 
Because of this probabilistic assignment Wolff clusters 
are subsets of the set of clusters formed by associating 
all adjacent spins with the same orientation. Clusters 
formed by considering spins with the same orientation, 
without this probabilistic restriction, percolate below Tc- 
It is the typical size of these latter 'conventional' clusters 
that is associated with defect densities, according to the 
Kibble-Zurek |p^ , pT| theory of defect formation. Our ob- 
servation of monopole percolation below Tc is compatible 
with this scenario. The interesting questions of testing 
this hypothesis and of determining the detailed relation- 
ship between Wolff clusters and topological excitations 
at criticality will be left for future work. 
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APPENDIX A: ENERGETICS OF LINE 
DEFECTS IN D-DIMENSIONS 

The free energy arguments of Sec. || can be general- 
ized to extended topological defects. In this appendix we 
consider the case of line defects. 

Let us look first at the case of one-dimensional defects 
N = D — 1. For a string- like object the core energy is 
proportional to its length Z, and the entropy can easily 
be calculated assuming that it behaves like a random 
walk. The number of different configurations for a closed 
random walk in a cubic Z?-dimensional lattice is given by: 



judge in this approximation, since the correction to the 
linear interaction is only logarithmic. In reality we have 



n = {2D)\atiI)-^''^ X r 



(Al) 



The first two terms count the number of possible closed 
random walk loops with I steps and the final factor takes 
into account the arbitrariness in the choice of initial point 
in the loop. We see that the entropy S — log(ri) grows 
linearly with I. The free energy is then (omitting dimcn- 
sionless constants): 

^^(0 ^al + Ei + T{j + 1) log(0 - \og{2D)Tl, (A2) 

where a is the string tension and Ei{l) the interaction en- 
ergy for the loop. In general Ej is a function not only of 
the string length but also of its detailed shape. We start 
by considering the case when the interaction energy can 
be neglected. Then we see immediately that the system 
will undergo a phase transition above a certain T, charac- 
terized by the proliferation of long strings This is because 
both the entropy and core energy have the same depen- 
dence on /, as was the case in the Kosterlitz-Thouless 
transition. 

This set of approximations gives a reasonable descrip- 
tion of the D — 2 Ising model, where there is no long 
range interaction between strings and probably also of 
strings in gauge theories in any dimension. In contrast 
to the case of monopoles, where the pair interaction was 
responsible for the confining phase at low-T, here it is the 
core energy of the string that keeps it at a finite length. 

The interaction energy is hard to estimate and depends 
in general on details of the underlying model as well as of 
particular configurations. We can make a rough estimate 
by assuming that the string remains a random walk and 
consequently that the distance between two points on 
the string is of order of l^^^. Then if two distant string 
segments interact as point-like defects in a 0{N~1), D—1 
theory, the loop energy will be given by: 



Ei{l)~ I X l'^^-''^/^ D>3 
Eiil)^nog{l), D = 3 



(A3) 



good evidence |32| that for 3D 0(2) there is a percola- 
tion transition in the string network at Tc at which long 
strings spanning the volume appear in the system. 

For higher dimensional defects the calculation of the 
entropy becomes much harder. For example surfaces can 
be topologically complicated, exhibiting holes and han- 
dles. The general trend seems to suggest that for fixed 
spatial dimension, the entropy of a typical defect config- 
uration decreases with N. Since the interaction energy 
increases with the number of field components, we expect 
that for any D, and above a certain N the free energy 
will always be dominated by the interaction component 
and proliferation of defects is prohibited at all finite tem- 
peratures. For Z? = 3, for example, strings proliferate in 
the high-T phase for N — 2 but no unbinding of charge 
pairs occurs for iV = 3. 



APPENDIX B: DEFINING WINDING CHARGE 
ON THE LATTICE FOR GLOBAL THEORIES 

We identify 0(3) monopoles in a 3-D cubic lattice by 
generalizing the well known 'geodesic rule' used in most 
cosmic-string lattice based simulations . 




Clearly for D > 3 the interaction term dominates over 
the entropy for all T indicating the absence of a prolifer- 
ating phase. The marginal case, for D = 3, is harder to 



FIG. 13. The monopole charge in a cubic cell is identified 
by projecting the field-vectors at every corner onto a unit 
sphere. Each triangle on the square faces cube is thus mapped 
into a a spherical triangle (the one with the smallest surface is 
chosen). The sum of the surface of all these triangles, divided 
by 4-K, is then taken as defining the monopole charge inside 
the cubic cell. 

We count the winding of the field- vectors around each 
unit cell in our grid by using a 'smallest area' assumption. 
To this end, we triangulate the faces of each lattice cube 
and then map the 0(3) field- vector at all corners onto the 
unit sphere (the presence of a monopole depends only on 
the orientation of the field, not on its norm). For each 
triangular element in the cube's surface this defines a 
solid angle on the unit sphere (see fig. ^3|). The sign of 
the solid angle is taken according to the handedness of 
the corners. Its value 8 can be calculated thanks to a 
formula that relates the area of the spherical triangle, 
defined by three vectors on a unit sphere, to the angles 



between the geodesic sides of the triangle: 
|e| = a + /? + 7-7r. 



(Bl) 



Summing the sohd- angles corresponding to all the 12 tri- 
angles in the cubes surface, we obtain ^ = 47r-n where 
n is an integer taken to be the charge of the monopole 
inside the lattice cube. \n\ has an upper bound of 5 but 
in practice we rarely observe charges larger that 2 (which 
we interpret as two coincident unit charges). 
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